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We evaluate the strangeness production from equilibrating 
and transversely expanding quark gluon plasma which may 
be created in the wake of relativistic heavy ion collisions. We 
consider boost invariant longitudinal and cylindrically sym- 
metric transverse expansion of a gluon dominated partonic 
plasma, which is in local thermal equilibrium. Initial condi- 
tions obtained from the self screened parton cascade model 
are used. We empirically find that the final extent of the par- 
tonic equilibration rises almost linearly with the square of the 
initial energy density. This along with the corresponding vari- 
ation with the number of participants may help us distinguish 
between various models of parton production. 
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I. INTRODUCTION 

The study of ultra-relativistic heavy ion collision has 
entered a new era with the successful commissioning of 
the Relativistic Heavy Ion Collider at Brookhaven. This 
provides an opportunity to verify the possible occurrence 
of a phase transition from hadronic matter to deconfined 
quark matter, where partons are the basic degrees of free- 
dom. At the (lower) SPS energies, an enhanced produc- 
tion of strangeness, considered to be one of the more 
robust signatures of quark-hadron phase transition has 
already been observed The initial temperatures 

likely to be attained at RHIC and the LHC are expected 
to be much larger. A natural question would now be: 
how quickly is the strangeness equilibrated, if at all? 

In heavy-ion collisions strangeness is produced abun- 
dantly through the partonic interactions if the tempera- 
ture T > 200 MeV, the mass threshold of this semi-heavy 
flavour. The extent of its equilibration would however 
depend upon such details as the thermal and chemical 
evolution of the partonic system and the life-time of the 
hot deconfined phase. It has recently been shown that 
the chemical equilibration of the light flavours and the 
gluons slows down due to the radial expansion ||,^. It 



should then be expected that the extent of strangeness- 
equilibration can also be affected if allowances are made 
for the likely radial expansion of the plasma. The present 
work attempts to get answers to this and related ques- 
tions. We limit our discussions to the production and 
evolution of strangeness during the deconfined phase, 
whose initial conditions are taken from the Self Screened 
Parton Cascade (SSPC) model which has formed the 
basis of a large number of related studies in recent times. 
An early work in this direction used the initial conditions 
obtained from the HIJING model and considered only a 
longitudinal expansion In the following we closely 
follow this treatment and extend it to include transverse 
expansion as well. 

Our paper has been organized as follows: Section II 
describes briefly the basic equations of the hydrodynamic 
and chemical evolution of the partonic gas through the 
partonic reactions in a (1+1) dimensional longitudinal 
expansion and a (3-fl) dimensional transverse expansion. 
A brief summary is given in section III. 



II. HYDRODYNAMIC EXPANSION AND 
CHEMICAL EVOLUTION 

A. Master Equations 

We start with the assumption that the early 
(semi)hard collisions among partons produce a thermal- 
ized partonic plasma. The high px partons produced 
early in the collision, then provide a colour screened en- 
vironment for the production of partons having low pT 
and the high density of the partons launches the Landau 
Pomcranchuk Migdal (LPM) suppression mechanism to 
eliminate the coUinear singularity in parton fragmenta- 
tion. This leads to the the so-called self screened parton 
cascade model [0| and can be used to provide plausible 
initial conditions for the system. The subsequent chemi- 
cal equilibration is then attained through reactions of the 
type gg ^ qq and gg ^ ggg. 

The evolution of the system is now controlled by the 
equation for conservation of energy and momentum of an 
ideal fluid: 
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, T'"' = (e + P)u^'u'' + Pg 



(1) 



where e is the energy density and P is the pressure mea- 
sured in the rest frame of the fluid. The four-velocity 
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vector u'^ of the fluid satisfies the constraint = —1. 

We assume that the distribution functions for partons 
can be scaled through equilibrium distributions as 



(2) 



where fj{Ej) = (e^^^ T 1)"^ is the BE (FD) distribution 
for gluons (quarks), and Xj (j — g, u, d, s) are the 



nonequilibrium fugacities, Ej — yPj+'^li ^^'^ ™j is 

the mass of the parton. 

Now one can write the number density, energy density 
and pressure for a_partially equilibrated multi-component 
partonic plasma 
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for i — u, d, s. We also have the colour and spin de- 
generacy di — 3 X 2 and Xi = rrii/T, where is the 
mass of the quark and K^s are modified Bessel func- 
tions. We take strange quark mass as 150 MeV. For 
massless quarks these expressions simplify considerably 
and we have 6i(0) = 2 • 9C(3)/27r2, 62(0) = 2 • 7n'^/A0 
and 63(0) = 1. We further assume that = A;^, which 
should be valid for negligible net-baryonic density. This 
should be a reasonable approximation at RHIC and LHC 
energies. The speed of sound (cg) can be obtained from 



dP 
de ' 



(5) 



We found it to be close to 1/ \/3, as we confine ourselves 
to T > 200 MeV. We must add that several lattice QCD 
evaluations suggest M that A = e — 3p > so that 



< 1/3. For such a situation the cooling of the plasma 
would be slower [^0| and thus a much larger time would 
be available for equilibration. 

We solve the hydrodynamic equations ([l]) with the as- 
sumption that the system undergoes a boost invariant 
longitudinal expansion along the 2-axis and a cylindri- 
cally symmetric transverse expansion pl l . It is then suf- 
ficient to solve the problem for z — 0. 

The chemical equilibration of the species j is governed 
by the master equation 



(6) 



where Rj are the rates which propel the system to- 
wards chemical equilibration. The system would be in 



chemical equilibrium when A, = 1 so that Rj{x) 



0. 



As mentioned earlier, the dominant chemical reactions 
through which the chemical equilibration proceed p^ ] 
are gg <-^- ggg and gg <-^- ii. Radiative processes involv- 
ing quarks have substantially smaller cross sections in 
perturbative QCD, and quarks are less abundant than 
gluons in the initial phase of the chemical evolution of 
the parton gas. Other elastic scattering processes ensure 
maintenance of thermal equilibrium. 

Under these assumptions the master equations for dif- 
ferent species become [p2| 



^^{ngu^') = (i?2^3 - i?3^2) - (^9-^ - ^^^9) 

i 



(7) 



in an obvious notation. 

The gain and loss term for the gluon fusion process, 
gg ii can be written as 



Rg^i ^ Ri^g — ( X^ — A,;A,) — 



d^pi 
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(8) 

where we have used the unitary relation \M.gg^{i\'^ = 
\M.ii^gg^ ■ The above integral can be written as free 
space cross section for fusion process folded with the dis- 
tributions for initial particles as 

igg^a J [(Jgg^UVl2] fg{pi)fg{p2) , (9) 

with the cross section given by 

1 f d^P3 f d^P4 



da 
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where vi2 = |vi — V2I, the relative velocity between the 
initial particles. Following Ref. ||l^ the Eq.(^) can be 
written 
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SO that the rates given in Eq.(g|) become 
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The net rate for the process gg <-> ggg can be written as 

1 f d^pi f d^p2 
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Similarly the integral in Eq.dl^) can also be written in 
factorised form [O 
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= Rang (1 - Xg) . (14) 
The density weighted rates in Eqs.( p^ , p^ ) are defined as 
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where the thermally averaged and velocity weighted cross 
section are 
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where superscript i in R2 and (T2 denotes to the fusion 
process for a given flavour. We give the details of these 
calculations in the next sub-section. In calculating the 
rates, R\ and i?3, we have also considered the tempera- 
ture dependence of the strong coupling constant as 
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with the cut-off Ao = 300 MeV. The scale, Q = 27rT, 
is derived by comparing the inverse of bare gluon prop- 
agator in the imaginary time formalism with the static 
longitudinal propagator p^ . 
Eq.(^ can now be simplified 
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= Rl 
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where Vr is the transverse velocity and 7 = 1/ ^ 1 — f ^ . 

If we assume the system to undergo a purely longitu- 
dinal boost-invariant expansion, Eq. (^ reduces to the 
well-known relation llSl 
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d^ 
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where r is the proper time. Using Eqs. (||) in Eq. (|T9), 
one can obtain the ultrarelativistic equation of motionlq] 

as 
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which reduces to 'D^{Q) = 63(0) = 1 for massless quarks 
(antiquarks). 

Now the master equations (^ can be written as 
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where 



Vnixi) = 
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(24) 

and I?„(0) = 1, for massless quarks. We can easily see 
that Eqs.(|l8|) reduce to Eqs.( |2^ , p^ ) if there is no trans- 
verse expansion so that the radial velocity, = 0. Now 
one can study evolution of the multi-component partonic 
plasma in terms of the parton fugacities by solving the 
Eq. (H^) for longitudinally and Eq. (|l^) for transversely 
expanding plasma with the parton chemical equilibration 
rates. 



B. Partons Equilibration Rates: 

In this subsection we briefly recall |^,|l3| the evaluation 
of equilibration rates for the gluon fusion process {gg — )■ 
a) and gluon multiplication process {gg — *■ ggg). 
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1. Gluon Fusion 

The differential cross section given in Eq.(|I^) can be 
written as 



d(j„ 



dt 



f 67r^(s, m\, m| 



(25) 



where ^ = s~ — 2s (m^ + ni^) + {m\ — rn^) , which re- 
duces to ^ = s^, for massless incident particles. The 
matrix element for the process, gg —^ ii, can be obtained 
in terms of Mandelstam variable from Ref. Il6|l as 
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= TT a. 
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in which Mi is the current quark mass. However, for 
massless quarks {q — u and d) the above Eq.(|2^) reduces 
to 



|2 = 16^2^? 



3 ut 
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This together with Eq.(p5|) diverges logarithmically as 
u, i — > 0. We assume that this logarithmic divergence 
for massless quarks can be regularized by assigning them 
mass given by the thermal mass p2[ 



'■u,d 



Ag + i (A„ + Arf) 



47r ^2 
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(28) 



For s-quark no such approximation is necessary and we 
use Eq.(|6|) with m,. = 150 MeV. 

Now integrating the matrix element in Eq. ( p5[ ) over the 
variable t, between the limits 



Now the thermally averaged, velocity weighted cross 
section is defined as 

= (o- jjui2) = 77— — — — — — . (31) 

J d-^Pld'^P2fg{Pl).fg{P2) 

Using vi2 = ^-^/^(s, 0, 0)/2pip2 and the thermal average 
of s for a pair of gluons, (s) = IST^, we get 



9 /C(2) 
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Combining Eqs.([l5|) and (^2|) we finally obtain the gluon 
fusion rate. 



2. Gluon Multiplication 

Following Ref. [|l2| we also estimate the for gluon 
multiplication process from the differential cross-section 
given by O 



dalf dug 



^ cosh 2/ 
dq\ d^k^ dy dq\ d'^k± dy \ k± 
X 6* (\/s — fcj^coshy) , 



(33) 



where the first step function includes the approximate 
LPM suppression of the induced gluon and the second 
step function accounts for energy conservation. Here k± 
denotes the transverse momentum, y is the rapidity of 
the radiated gluon and gj_ corresponds the momentum 
transfer in the elastic collisions. The infrared divergence 
associated with QCD radiation is regularized by the LPM 
effect. However, Eq. (|33| ) is still having infrared singu- 
larities in both scattering cross sections and radiation 
amplitudes associated with the gluon propagator. One 
can, approximately, control all these singularities using 
the Debye screening mass [O] 



Qq2 foa 

= ^ / kfg{k)dk = AirasT'^Xg 
Jo 
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Now the regularized gluon density distribution induced 
by a single scattering is 
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one obtains the total cross-section 
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where x = ^1 — AMf / s. 



M^\ 1 + Y 
J 1 - X 



(29) 



(30) 



and the regularized small angle gg scattering cross section 
is 



do-ff _ 9 27ra2 



(36) 



The mean free path for the elastic scattering is obtained 
as HI 



Aj = —aiasl - 



8 l + 8na,\g/9 



(37) 
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Integrating the (/) part analytically we get the gluon mul- 
tiplication rate as 

32 / 8 \^ 

R3 = ^^(^sTXg h + -aiasXgj J{\g) , (38) 



where 



1{\g) = dx 



dz- 



cosh""^ (a/x) 
x\J[x + (1 + z)XD^ — AXZXD 

1 cosh"""^ (v^) 
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(39) 



D 



with Xd = Mzj-^/ = ^I'jj/s. 

In the following subsection we discuss the chemical evo- 
lution of different parton species. 



C. Evolution of the Multi-component Parton Plasma 

As mentioned earlier we shall be using the initial con- 
ditions obtained from the self-screened parton cascade 
model . Even though they are well-known and used in 
a large number of studies we reproduce them in Table-I 
for the sake of completeness and an easy reference. We 
have chosen an initial fugacity for the strange quarks as 
half of that for the light quarks as in Ref. Q . This is con- 
sistent with the assumption, which is often made, that 
the number of flavours is ~ 2.5 if the mass of s quark 
is taken as zero [|2[. We show our results for the lon- 
gitudinal expansion for RHIC (upper panel) and LHC 
(lower panel) energies in Figure 1. As the additional 
parton production consumes energy, the temperature of 
the partonic plasma drops considerably faster than the 
ideal Bjorken's scaling (T = To{To/Ty^-^ , Tq and tq, re- 
spectively, are initial temperature and time of the parton 
gas) represented by the dashed line. We see that like the 
light quarks and gluon, the production of strange quarks 
continues till late in the evolution. We further note that 
the extent of equilibration for the strange quarks in com- 
parison to that for the light quarks (Xs/^u.d) rises rapidly 
and once the temperature falls below about 300 MeV 
('^ 2tos ) it gets more or less frozen by this time. Thus 
we conclude that the equilibration of strangeness pro- 
duction may imply the existence of quark matter at a 
temperature of more than about 300 MeV for a time of 
the order of a few fm/c. The other aspects of varia- 
tion of Ag, Xu.d and T have already been discussed by 
several authors We do note that the plasma is 

not fully equilibrated chemically at either RHIC or LHC 
energies by the time the temperature drops below 200 
MeV. We do expect that additional quark production 



may occur due to gluon fragmentations during hadroni- 
sation |^ 19,2^] leading to a chemically equilibrated hot 
hadronic matter. 

Before giving our results for transverse expansion we 
need to specify the profile of the initial energy den- 
sity and the fugacity of the partonic system. Following 
Ref. [|j2l]j2| we take the initial energy density as 



1-^ 

i?2 



1/2 



e(i?-r), 



(40) 



where R is the transverse dimension of the system, r is 
the radial distance, and Eq is the "average" initial energy 
density (see Tab. I). The profile plays an important role 
in defining the boundary of the hot and dense deconfined 
matter and affects the transverse expansion through the 
introduction of pressure gradients. We have further taken 



Aj(r, To) oc £(r, Tq) as before [|6|, pi| , p2| . Any other varia- 
tion will require an additional parameter. We give our 
results for the radial variation of Ag, Xu^d and Xs/Xu,d in 
Figs. 2 and 3 for RHIC and LHC energies, respectively, 
for various times along the constant energy density con- 
tours with T = Ntq. Here N is defined Js],^ through 
£(r,r) =e(r = 0,To)/iV4/3. 

We sec that the fugacities attain their highest values 
near r = and rise rapidly first and only slowly later 
in time. We also see a result unique to chemical evolu- 
tion with transverse expansion that the fugacities may 
even start reducing towards the end of the QGP phase 
when the radial velocity (gradient) becomes very large 
(see Ref. for explanation). 

We can use our results for the radial variation of the 
energy density and the fugacities obtained here to esti- 
mate the extent of partonic equilibration as a function 
of the initial energy density. Comparing results for this 
from figures 4 and 5 we note that once the energy den- 
sity is beyond about 20-40 GeV, the final fugacities for 
all the partons increase almost linearly with the square 
of the energy density. This is a very interesting result. 

Let us identify the "local" energy density at the ra- 
dius r with the average energy density attained in a 
non-central collision or with a central collision involving 
lighter nuclei. Recall that the partonic models suggest 
that the initial temperature attained in such collisions 
through the production of minijets varies as A^^^ p3| , p4t . 
This should then imply that the extent of equilibration 
of strangeness produced would rise as N^[^^, if we can 
identify A with 7Vpart/2, where A'part is the number of 
participants. In actual cases the variation may be some- 
what modified from this naive expectation due to the 
considerations of shadowing and jet quenching [ p5| . 

If on the other hand the energy densities etc. are de- 
cided by the considerations of parton saturation |2^ then 
the initial temperature would vary as ~ ^0.126 g^j^^ ^j^g 
energy density as ^ A'^-^. The results of Figs. 4 and 
5 would then imply that the extent of strangeness equi- 
libration should increase as the number of participants. 
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^part (see also Ref. for a slight reformulation of the 
parton saturation model). 

It is also interesting to note that a recent analysis |2^] 
of centrality dependence of the extent of strangeness equi- 
libration at CERN SPS energies in Pb+Pb collisions 
gives a linear increase for this with the number of partic- 
ipants. 



on the initial conditions and also on the evolution mech- 
anism. 



Wong has studied |30| the chemical equilibration of 
plasma essentially in a formalism similar to that used in 
the present work, but only with longitudinal expansions. 
His values for as are also much larger. 



D. Comparison with other Works: 

Rafelski and coworkers have studied the strangeness 
equilibration for more than two decades in various de- 
tails. We shall refer to only a recent work by these au- 
thors In Ref. 1^ the following assumptions have 
been made which are at variance with our works: 

1. The light quarks and giuons arc assumed to have 



attained chemical equilibrium (Ag = = 
when the evolution of strangeness begins. 



A. 



1 



2. A Boltzmann approximation is used to describe the 
phase space distribution of s-quarks. 

3. The flavour changing process, qq ^ ss, has been 
included. We have neglected them assuming that 
the gluon fusion process gg ss dominates in a 
chemically undersaturated plasma. 

Now in the Boltzmann limit the evolution of the s- 
quark fugacity can be written from (|2^) as 



1 



Xs Ki{Xs) 



3 K2{xs) 
This can be rewritten as 



r = ^Afi-^)- (41) 



Tfis 



It 



Xs_ Ki{xs 
T'^'K^ixs 



3^ 



T 



(42) 



which corresponds to Eq.(ll) of Ref. |2£|. We see that 
the 3''^ and 4'^ terms in the left hand side here are absent 
in Eq.(ll) of Ref. |^^, due to the above discussion. One 
can thus see that the extent of strangeness equilibration 
attained in the work of Ref. ||2^] is much larger than in 
the present work. 

As a further check, we have repeated our calculations, 
assuming Ag = A„ = A^; = 1 with the same energy density 
as before, (see Table I) and further taking As(to) = 0.2 as 
in the work of Ref. |29|] . We found that now our parame- 
ter Xs rises to about 0.4 at RHIC energies and up to 0.72 
at LHC energies, for the case of longitudinal expansion 
(see Fig. 6). When the transverse expansion is allowed 
these numbers reduce to 0.66 at LHC energies and 0.38 
at RHIC energies, respectively. We thus realize that the 
extent of equilibration of strangeness depends sensitively 



III. SUMMARY 

We have studied the evolution and production of 
strangeness through the partonic interactions in a chem- 
ically equilibrating and expanding multi-component par- 
tonic gas with the initial conditions obtained from SSPC 
model. We find that most of the strange quarks are pro- 
duced when the temperature is still more than about 300 
MeV (~ 2ms). Thus a chemically equilibrated plasma 
is expected to imply the existence of QGP phase for a 
duration of several fm/c. We also find approximately 
that the extent of strangeness equilibration rises linearly 
with the square of the initial energy density within our 
approach. This may help us to obtain the scaling of the 
initial energy density with the number of participants and 
distinguish between the minijet and the partonic satura- 
tion models of parton production. It is rather interesting 
that the charged particle multiplicity in Au -I- Au col- 
lision at y/SNN = 130 GeV measured by the PHENIX 
collaboration [ pl| shows a behaviour which is a superpo- 
sition of two terms, a linear increase with the number 
of participants and a linear increase with the number of 
collisions (which varies as N^^^^). 
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TABLE I. Initial conditions for the hydrodynamical ex- 
pansion phase in central collision of two gold nuclei at BNL 
RHIC and CERN LHC energies, respectively, from SSPC § 
given in first and second row. Initial conditionis given in third 
and fourth row are obtained with the same energy densities 
as before along with A„ = A„ = Ag = 1. 



Energy 


To 

(fm/c) 


To 
(GeV) 


Af 


x(o) 

u[a.) 


Ai"' 


(GeV/fm^) 


RHIC 
LHC 


0.25 
0.25 


0.67 
1.02 


0.34 
0.43 


0.068 
0.086 


0.034 
0.043 


61.4 

425 


RHIC 
LHC 


0.25 
0.25 


0.44 
0.72 


1.0 
1.0 


1.0 
1.0 


0.2 
0.2 


61.4 

425 
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2 4 6 8 10 12 14 

T [ fm/c ] 

FIG. 1. Variation of temperature, coupling constant, gluon 
and quark fugacitics with proper time for (1+1) dimen- 
sional hydrodynamic expansion with SSPC initial condition 
for RHIC (upper panel) and LHC (lower panel) energies. 
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FIG. 2. Variation of gluon (upper panel), massless 

quark (middle panel) fugacities and ratio of strange to 
nonstrange quark fugacity (lower panel) with the trans- 
verse radius for RHIC energy at different times r = Nm, 
along the constant energy density contours defined by 
e(r,r) =£(r = 0,ro)/iV*''^ 
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FIG. 4. Variation of the final fugacities with the initial en- 
ergy density for RHIC energy. Note the scale of the x-axis 
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FIG. 5. Same as Fig. 4 for LHC energy. Note the scale of 
the X-axis. 



FIG. 3. Same as Fig. 2 for LHC energy 
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FIG. 6. Sensitivity of the strangeness equilibration to the 
initial conditions for longitudinal expansion. The solid curves 
give the results corresponding to Fig. 1 given earlier, while the 
dashed curves correspond to initial condition which have the 
same energy densities as before, but have = A„ = Ag = 1 
and As = 0.2, as given in Table I. 
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